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ABSTRACT 

The seasonal variations in temperature, salinity, water 
masses and geostrophlc currents at the surface and 200 m are 
discussed for the region off the Central California coast 
near Monterey. The results were obtained from data collected 
on six cruises covering a span of two years from 1970 to 
1972. 

It is shown that the surface curren1:s are subject to 
marked seasonal variability with the mean flow characterized 
by eddies and filaments. The currents at 200 m are similar 
to those found at the surface. The bottom topography 
influences the direction of flow inside the 1000 fathom 
curve . 

Seasonal variations in temperature and salinity generally 
conformed to Skogsberg’s (1936) annual cycle composed of 
three distinct hydrographic seasons. 

Only subtle changes in water masses from season to season 
were noted. Greater than 60 per cent Equatorial Pacific 
v;ater is found throughout the year below 800 m. The propor- 
tion of Equatorial Pacific water increases with depth and 
towards shore . 
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I . INTRODUCTION 



Kln'dyushev (1970) stated that the California region is 
one of the most intensely studied regions of the Pacific 
Ocean. This statement may be factual, but only with certain 
reservations. A review of the available literature reveals 
that the region from the Gulf of Alaska south to Cape 
Mendocino (4l°N) has held the Interest of researchers for 
many years. This may be explained by: (1) the uniqueness 

and complexity of the oceanic processes that occur in this 
region; (2) the high concentration of commercial fisheries 
and other important commercial enterprises; and (3) the 
fact that several oceanographic institutions such as Oregon 
State University and the University of V/ashlngton are 
physically located adjacent to the area. The region from 
Point Conception (3^°N) south to Cabo San Lucas on the 
southern tip of Baja California has received a great deal 
of attention by Scripps Institute of Oceanography and 
various less well known institutions for the same reasons 
cited above. These constitute the areas of Intense research 
along the west coast of the United States. Thus there is 
a region along the Central California coast between Cape 
Mendocino and Point Conception that has received relatively 
little attention in the literature. 

The purpose of the research reported in this study is 
to help close the "gap" noted above and to document an 
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analysis of the geostrophlc circulation patterns off the 
Central California coast near Monterey. The hydrography of 
the area is examined and an attempt is made to determine the 
distribution of the various water masses that are present 
in this region. An additional objective was to examine the 
validity of the seasonal concept for this region first 
postulated by Skogsberg (1936), wherein the annual cycle is 
divided into three distinct hydrographic seasons: the 

Davidson Current, Upwelling, and Oceanic Seasons. 

The region off the Central California' coast is charac- 
terized by complex flow patterns that are subject to diurnal, 
semidiurnal and seasonal variations (Reid and Swart zlose, 
1963 ). The offshore California Current, the inshore Davidson 
Current, seasonal upwelling and the bottom topography are 
just some of the factors which play an important role in the 
formation of the currents in this region. In addition, the 
area lies in a zone in which both subarctic and subtropical 
water masses are present. This is known as the "transition 
zone," a term first introduced by Sverdrup, Johnson, and 
Fleming (19^2). As noted by Kin'dyushev (1970), considera- 
tion of current patterns in an area is essential for an 
understanding of the origins and interactions of the various 
waters observed in that region. 

Plow in the surface layers of the deep ocean is influenced 
by several forces: the most important are the horizontal 

pressure gradient and viscous forces within the fluid, the 
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action of wind stress, tidal forces, and the Coriolis force. 
At present our incomplete knowledge of ocean dynamics and 
the lack of suitable available data, such as direct current 
measurements, restrict this study to the indirect deter- 
mination of surface currents. Accordingly the "dynamic 
method" for computing the geostrophic current is utilized. 
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II. REVIEW OF THE LITERATURE 



A. GENERAL CURRENT PATTERN 

The currents of the eastern north Pacific Ocean are 
part of the large anticyclonic gyre centered near the 
Hawaiian Islands. The eastern portion of this gyre — the 
California Current — represents the continuation of the 
Aleutian Current of the North Pacific (Sverdrup, Johnson 
and Fleming, 19^2). As can be seen in Figure 1 it develops 
as an Independent current in the region of divergence of 
the North Pacific Current (approximately 48®N) at the North 
American coast. It flows southeast along the continent at 
0.2-0. 6 knots in a band approximately 400 miles wide 
(Kin* dyushev, 1970). Due to the large expanse of the 
current, no high velocities are encountered except within 
local eddies (Sverdrup, et al., 1942). In the vicinity of 
latitude 32°N the California Current turns toward the coast 
and then to the north in a sharp meander before turning 
south again. Pavlova (1966) reported that this stable 
rotation is probably due to features of bottom topography 
in the area and to the pattern of prevailing winds along 
the coast. Pavlova's results were based on analysis of 
geostrophlc currents. Since then the rotation of the water 
at 32°N was confirmed with parachute drogues by Reid, 

Swartzlose and Brovm (1963). Pavlova also noted that 
rotations could be traced on charts of geostrophlc circulation. 
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Fig. 1 . Mean surface currents in the North Pacific (from Namias, 1970). 
Solid curve indicates trajectory of water from A to B around Pacific gyre. 
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Generally speaking, south of 32®N the California Current 
represents a wide body of water which moves sluggishly 
toward the southeast and eventually turns westward at about 
24°N (Reid, Roden, and Wyllie, 1958). 

Characteristic of the California Current are significant 
seasonal variations (depending on changes in the wind field) 
associated with the seasonal displacements of the center of 
the Hawaiian anticyclone (Reid, et al., 1958). During the 
fall and winter months, when the prevailing northerly winds . 
weaken, a countercurrent, referred to by many authors as the 
California Countercurrent, or more commonly as the Davidson 
Current, develops in the surface layer between the California 
Current and the coast. At this time the Davidson Current is 
most intense and takes the form of a narrow band within this 
region. In spring and early summer north-northwest winds 
intensify and prevail over the entire coastal region 
(Sverdrup, et al., 19^2). During these months the prevailing 
winds give rise to intense upwelling and the Davidson Current 
gradually breaks up into separate cyclonic and anticyclonic 
eddies (Kin ’ dyushev, 1970; Sverdrup, et al., 19^2). Pavlova 
(1966) reported that north of Point Conception, the Davidson 
Current reaches its maximum development at depth (200-250 m) 
in the summer and autumn. It reaches its maximum development 
at the surface from October through April, i.e., late autumn 
to early spring. - Furthermore, in August the Davidson Current 
is scarcely noticeable at the surface despite active develop- 
ment at depth. In December the core of maximum velocity 



13 



emerges at the surface and In the late spring it almost 
completely disappears (Reid, I96O; Reid, et al., 1958, 
Pavlova, 1966). 

Based on data collected within the confines of Monterey 
Bay, Skogsberg (1936) divided the year into three periods: 
the Davidson Current season, lasting from the middle of 
November to the middle of February; the upwelllng season, 
between the middle of February and the end of July; and the 
oceanic season, from the end of July to the middle of 
November. He characterized the Davidson Current season as 
one in which there is a relatively thick upper homogeneous 
layer where the temperature at 50 m is only slightly lower 
than at the surface. During the upwelling period the 
California Current reaches the coast and, under the influence 
of the prevailing northwesterly winds, upward vertical motion 
takes place near the coast bringing relatively cold sub- 
surface water toward the surface. As a result, the tempera- 
tures at given depths decrease when upwelling begins. As 
the upwelling gradually ceases toward the end of August, a 
sharp rise in temperature takes place both at the surface 
and at subsurface depths due to effects of heating and 
conduction and intrusion of offshore water. 

B. WATER MASSES 

A number of special studies of the water masses of the 
California region of the Pacific Ocean have been reported 
by Tlbby (19^1), Sverdrup and Fleming (19^1), Sverdrup, 
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et al., (19^2). Reid, Roden and Wyllie (1958) described 
the hydrography of the entire Northeast Pacific Ocean and 
reported that water entering the California Current system 
comes from four great water masses. The offshore waters 
of the northern part of the California Current are derived 
from the Subarctic Water Mass. As this water moves to the 
south (east of 1^7°W), it mixes with waters from the Central 
Water Mass which move in from the northwest and west. A 
major influx of southern water from the Equatorial Water 
Mass enters the system as a subsurface current from the 
south. This current is the previously discussed Davidson 
Current which tends to be near the coast. The fourth major 
water source is derived from the upv;elling processes at the 
coast . 

Sverdrup and Fleming (19^1) > Tibby (19^1) > and Sverdrup, 
et al, (19^2), simplified the classification problem by 
defining two extreme water masses: a "northern” and a 

"southern." These two water masses were named according 
to the regions from whence they originated and were called 
"Subarctic North Pacific" and "Equatorial Pacific," 
respectively . 

A characteristic of Equatorial Pacific water, as defined 
by Sverdrup and Fleming (19^1), v;as a Temperature-Salinity 
(T-S) curve on which salinity was relatively constant as 



The term, v/ater mass, is used in the manner defined 
by Sverdrup, Johnson, and Fleming, 19^2: 1^<3-1^6. 
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temperature decreased (see Figure 2, curves S.3 and B-III, 
31 ). A T-S curve for Subarctic water shows an Increase in 
salinity with decreasing temperature (Figure 2, curve C 131 ). 
They, as well as Tibby (19^1), constructed a chart defining 
percentage Equatorial water for given T-S pairs (Figure 3). 
Based on this definition they traced the northward intrusion 
of Equatorial water. As might be expected, the percentage 
of Equatorial water decreased in the direction of flow in 
the northward flowing current . 

The region of study off Monterey lies nearly in the 
center of the transition zone between the subarctic and the 
subtropical water masses. Roden (1971) placed the northern 
boundary of the transition zone near ^12°N, and the southern 
boundary at approximately 30°N. However, he, as well as 
most other investigators, concentrated on determining the 
boundaries of the transition zone rather than on conditions 
within the zone. 
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Figure 2. Temperature-salinity curves 
from selected stations (from Sverdrup 
and Johnson, 19^1). 




Figure 3. Graph shov/lng T-S curves 
defining Subarctic North Pacific and 
Equatorial v/ater, and curves for various 
percentages of Equatorial v/ater assuming 
mixing along surfaces of equal 
(from Tibby, 19^1) • 
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Ill . THE DATA 



A. AREA OF INVESTIGATION AND COLLECTION OF DATA 

The area of investigation is sho^m in Figure in 
which the outstanding features of the topography of the 
bottom have been indicated. The major topographical 
features of the area include the Monterey Canyon and 
Monterey Bay. The area of investigation is bounded on the 
north by latitude 37°N, on the west at approximately 124°W, 
and on the south by latitude 35°30'N. 

The data were obtained from six U.S. Naval Postgraduate 
School oceanography student training cruises. The station 
plans for the various cruises are shown in Figures 5-10. 

It should also be noted that the cruises were conducted 
during varying months, over a period of two years, thus 
making available data for the study of seasonal variability 
of properties and conditions. 

Water samples were obtained by means of standard Nansen 
reversing water bottles, and temperatures by means of 
reversing thermometers. On each cast both unprotected and 
protected thermometers were attached to tv;o or three Nansen 
bottles in order to accurately determine the depths of 
observation. All other bottles were provided with tv/o 
protected thermometers. 

The temperatures have been corrected in the ordinary 
manner. The salinities viere determined on board the 
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Fig. 4. Major geographical and, topographical features of the study area. 
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Pig. 5. Oceanographic stations occupied during Cruise 
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Fig. 6. Oceanographic statiorts occupied during Cruise II. 
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Fig. 7. Oceanographic stations occupied during Cruise III. 
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Fig. 8. Oceanographic stations occupied during Cruise IV. 
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Fig. 9. Oceanographic staltlons occupied during Cruise 
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Fig. 10. Oceanographic stations occupied during Cruise VI. 



research vessels using the Bis sett -Berman model #9006 
salinometer. The accuracy of the salinity determinations 
is generally on the order of i 0.01 %© . The data for all 
stations are contained in Appendixes B-G, which furthermore 
contain the computed values of Sigma-t. The values were 
computed by means of program "HYDRO" developed for use with 
the IBM 360 computer by the U.S. Navy Postgraduate School 
Department of Oceanography. Program HYDRO is presented in 
Appendix A. 
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IV. ANALYSIS 



A. PRESENTATION OF DATA 

Considerable time and effort were spent on a critical 
review of the temperature and salinity data. This was 
considered essential since the data had been collected 
mainly by inexperienced personnel. Temperature-salinity 
curves were prepared for each station and used for dis- 
covering possible conspicuous errors. In addition, these 
T-S curves were used to determine the percentage of 
southern water in the area during different seasons of the 
year. The distribution of southern water is presented 
graphically in vertical sections along transects approxi- 
mately normal to the coast. 

Vertical cross sections of salinity and temperature 
along latitudinal transects were prepared for the purpose 
of determining the distribution of properties in the area, 
and as an aid for detecting upwelling. In addition, charts 
showing the horizontal distribution of temperature and 
salinity at the surface and at 200 m were prepared for the 
same reasons cited above. 

Charts depicting geostrophlc flow at the surface; and at 
200 m were prepared using 600 m as the reference level. 

The dynamic method for computing geostrophlc flow and the 
selection of the reference level are discussed in the next 
section. It should be noted that the geostrophlc flov; 
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patterns are represented by lines of mean flow Instead of 
contours of dynamic height anomolles which are used by 
many authors. Small numerals Integral to the flow lines 
Indicate current speeds In cm/sec. 

B. THE GEOSTROPHIC CURRENT 

Geostrophlc computations, also often referred to as 
dynamic computations or the dynamic method, are discussed 
In any elementary physical oceanography text, for example, 
Neumann and Pierson (1966). Fomin (1964) discusses in 
great detail the application and the limitations of the 
dynamic method. What might be called the classic computa- 
tional scheme Is based on an equation originally derived 
from Bjerknes circulation theorem (1900) by Sandstrom and 
Helland-Hansen (1903) . However, the resulting equation Is 
commonly referred to In the literature as the Helland- 
Hansen equation and It has been used extensively In ocean- 
ography to compute the relative field of currents from the 
observed field of mass. The equation, written In terms of 
the horizontal pressure gradient of the geopotential between 
two levels P^ and ?2 Is, 
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where 



f = 2ti) sin (f) is the Coriolis parameter (1/sec) 

(I) = the angular velocity of the earth (radians/sec) ’ 

(j) = the geographic latitude (degrees) 

a = specific volume C(cm^/gm), subscripts A and B refer 
to stations A and B separated by distance AX] 

(V 2 -V^)= the horizontal velocity at relative to the 
level P 2 (cm/sec) 

AX = the distance between stations (kilometers) 

P = pressure (decibars) 

The geopotential difference between two levels in the 
ocean is determined by the integrals in the right-hand term 
of Equation 1. The integrals are numerically evaluated from 
vertical density distributions, determined from the measured 
vertical temperature and salinity structure and the equation 
of state. The measured temperature and salinity structure 
at two stations can be used in this way to compute the mean 
horizontal gradient in geopotential along an isobaric 
surface. Then Equation 1 can be evaluated to give the mean 
geostrophic flow normal to the horizontal line joining the 
two stations. 

From the observed mass distribution only the relative 
geopotential with respect to the sea surface and hence the 
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relative velocity field can be derived. In order to arrive 

at the absolute current, either the absolute geopotentlal 

field or the absolute velocity must be known at least at 

one level between two adjacent stations (Neumann and 

Pierson, 1966). It Is common practice to satisfy this 

requirement by assuming a depth where the currents become 

zero or by Inferring It from considerations of continuity. 

This depth Is commonly referred to as the level of no 

motion or the reference level. 

What constitutes the proper reference level Is open to 

question (Fomin, 1964). Many reference levels have been 

used to compute surface currents (Neumann and Pierson, 1966). 

Reid ( 1958 ) used a reference level of 500 db In the 

California Current region. Sverdrup and Fleming (1941) also 

used a 500 db reference level for dynamic computations In a 

section off Southern California. In a study of the 

California Current regime north of San Diego, Reid, Swart zlose 

and Brown (1963) showed that geostrophlc currents computed 

with respect to a 500 db reference level compared favorably 

with surface currents Independently measured by drogues. 

Stommel (I 965 ) summarized the problem of selecting a 

reference level (page 20): 

Thus the choice of the reference level for 
geostrophlc calculations becomes mostly a matter 
of taste, and v/e should admit that that Is 
ultimately Intolerable. The determination of the 
level of no motion is not a matter for debate, but 
for direct measurement . 
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There are simply not sufficient direct current measure- 
ments in the oceans to select suitable reference levels. 
Adding to the complexity of the problem is the fact that 
even if there were direct measurements in the ocean, only 
those in the particular flow regime under study would be 
useful. Fomin (196^) shows that in cases where the 
horizontal velocity approaches zero at some depth and remains 
nearly so with increasing depth, it may be better to select 
the shallower reference level of low velocity rather than a 
substantially deeper level, thus reducing computational 
errors. He point out that in some cases the selection of a 
deep reference level may completely distort the computed 
currents. For this study a reference level of 600 m was 
selected which permitted a substantially greater number of 
stations to be used. Stommel ( 1965 ) states (page 19 ) » 
"Fortunately, the choice of reference level has less effect 
upon velocity in the near surface layer than on those of 
deep water." 

C. LIMITATIONS OF THE DYNAMIC MODEL 

Use of Helland-Hansens ' formula to infer absolute 
currents suffers from several limitations (Fomin, 196^). 

Use of the dynamic method for this is based on the 
follov/ing assumptions: 

( 1 ) The flov; is geostrophic (frictionless and 
non-accelerated) . 
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(2) The horizontal pressure gradient force and 
the horizontal velocity become negligible at a moderate 
depth below the sea surface. In addition, either the depth 
where the horizontal velocity becomes negligible must be 
known, or the horizontal velocity remains negligible with 
increasing depth. 

Newmann and Pierson (1966) warn against the uncritical 

use of the dynamic method. They state (page 182), 

Nothing more can be expected of this method than 
its simple assumptions permit. This, however, 
is true for all theoretical approaches, and 
these critical remarks neither inflate nor 
depreciate the value of simple analytical attempts 
such as the "classical method of dynamic 
computations," so long as their limitations or 
restrictions are kept well in mind. 

However, experience has shown the dynamic method to be 

highly useful and "reliable" in deep oceans. One outstanding 

example of the computation of an absolute current profile, 

using the dynamic method is a study of the currents in the 

Straits of Florida by Wiist (192^), (as cited by .Sverdrup, 

Johnson and Fleming, 19^2) . 
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V. RESULTS 



A. REPRESENTATION OF CURRENTS 

In the treatment of the data the usual convention of 
first describing the temperature and salinity distribution 
will not be followed. Instead, the results of the dynamic 
computations are presented first, since a discussion of 
the water masses will be simplified if the character of the 
flow has first been described. In addition, the results 
are presented in a logical order according to season of year 
as follows : 

Winter (or Davidson Current Season) 

Cruise I 8-14 November 1970 

Cruise II 19-23 January 1972 

Cruise III — 4-17 February 1970 
Cruise IV 27 February-4 March 1971 

Early Summer (or Upwelling Season) 

Cruise V 27 June-2 July 1971 

Late Summer (or end of Upwelling Season, beginning 
of Oceanic Season) 

Cruise VI— 8-14 August 1970 

It is recognized that there are uncertainties which 
arise v/hen attempting to drav; conclusions about the 
circulation computed from the dynamic topography. Of major 
importance in this regard is the influence of internal 
vmves. Defant (1961) points out the fact that internal 
waves distort the details of the mass field by periodically 
displacing the water masses in a water column. As a result. 
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these periodic variations in ocean structure will show in 
the dynamic evaluations. However, it cannot be doubted 
that the currents derived from dynamic computations 
represent the general direction of the flow. For example, 
a broad view of the seasonal variation is shown by the 
mass transport (see Table 1, pageil3 ), which agrees with 
gross features of the flow system shown for the corresponding 
cruises by the surface current charts. Prom the dynamic 
topography an approximately correct picture of the flow can 
be obtained, and an attempt can be made to verify this by 
a close study of the characteristics of the different water 
masses. It should be noted that construction of current 

I 

charts is, in general, subjective. Due to the complexities 
involved in this current regime, they may be subject to 
more individual interpretation than is ordinarily required. 

As stated previously, the currents are represented 
schematically by lines of mean flow Instead of by contours 
of dynamic height anomolies. Small numerals integral to 
the lines of flow indicate the current speed in cm/sec. 

B. SURFACE CURRENTS 

Surface current patterns during winter (Cruises I, II, 

III and IV) are shown in Figures 11, 12, 13 and 1^. The 
surface currents were similar in one respect during all 
four cruises. Near the coast, south of Monterey Bay, was 
a well defined current to the northv/est which ran approxi- 
mately parallel to the coastline. This was expected since 
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Fig. 13. Cruise III. Surface currents (relative to the 6OO db reference level) 
represented by solid lines with arrows indicating direction of flow. 

Current speeds are in cm/sec. 
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the winter season is normally the period when the Davidson 
Current reaches its maximum intensity at the surface. 

Of greater interest was the variability of the flow 
throughout the winter season. The November 1970 cruise 
results, which depict conditions at the start of the winter 
season, show that the California Current is still evident 
relatively close to the coast near the 1000 fathom curve 
at an average speed of approximately 10 cm/sec. The near- 
shore poleward current was somewhat faster, averaging nearly 
20 cm/sec. and appeared to decelerate along the direction of 
flow. Results of the January 1972 and February 1971 cruises 
are remarkably similar, showing a well defined poleward 
offshore current and two distinct eddies; one near Monterey 
Bay, and the other just north of Pt . Sur. The northern 
eddies were cyclonic on both cruises. However, the southern 
eddy in January 1972 was anticyclonic , whereas the southern 
eddy in February 1971 was cyclonic. The offshore poleward 
floxif was considerably higher in January 1972 (Cruise II), 
averaging nearly 20 cm/sec, whereas a relatively weak 
poleward current averaging less than 10 cm/sec was noted 
in February 1971 (Cruise IV). The chart for the February 
1970 cruise shov;s several interesting details not evident 
on the other cruises. Overall weak polev/ard flow is 
Indicated, even as far , out as 110 nm from the coast. 

However, approximately 50 nm offshore a band of equatorward 
flow interlaces betv.’een separate bands of northerly flow. 
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This finger of current to the south appeared to accelerate 
in the direction of flow and was separated from the poleward 
flow on its western side by an anticyclonic eddy. Station 
spacing was not sufficiently dense to conclude whether small 
scale eddies were present on the eastern side of the 
equatorward current, or if simple lateral shear separated 
the oppositely directed currents. 

Of particular interest on all four cruises was the 
observation that the bathymetry appeared to exert a definite 
Influence on the current direction, especially coastward of 
the 1000 fathom curve. 

Figure 15 shows the surface currents during early summer 
(Cruise V) for which overall weak, poleward flow is evident. 
However, the area is characterized by small scale cyclonic 
eddies (on the order of 10 miles) .indicating the general 
degradation of the intensity of the Davidson Current at the 
surface. This appears to agree well with Kin ’dyushev’ s 
(1970) observations. Also of interest is the finger of 
relatively strong equatorward flow between the nearshore 
eddies and the coast. Water mass analysis shows that this 
water is predominantly of Subarctic origin. This corrobo- 
rates Skogsberg’s (1936) observations that the California 
Current reaches the coast during the upwelllng season. 

Figure l6 shows the surface currents found on Cruise VI 
conducted during the late summer. The California Current 
was clearly evident at the stations farthest offshore. 
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Fig. 15. Cruise V. Surface currents (relative to the 600 db reference level) 
represented by solid lines with arrows indicating direction of flow. 

Current speeds are in cm/sec. 
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The Davidson Current was still dominant in the nearshore 
region and appeared to increase in speed along the direction 
of flow. Two large antlcyclonic eddies were indicated; one 
near the mouth of Monterey Bay, and the other farther off- 
shore south of Pt. Sur. Again, as with all previous cruises, 
bathymetry appeared to be an Influential mechanism inside 
the 1000 fathom curve . 

Table 1 is a compilation of volume transports between 
selected stations for Feb. 1970, June-July 1971, and Aug. 

1970 cruises. The volume transport may be interpreted 
as Indicating the predominant direction of flow between the 
stations; negative values indicating southward flow, and 
positive values indicating poleward flow. In order to form 
the best possible basis for comparison, matching station 
pairs of similar location were selected for each cruise. 



CRUISE III 


CRUISE V 


CRUISE VI 


FEB. 1970 


JUNE-JULY 1971 


AUG. 1970 


Stations Q and S 


Stations III-5 & III-9 


Stations 19 and 23 


AX = 6l km 


AX = 52 km 


AX = 66 km 


- 0.280 


+2.160 


+ 2.564 



6 ^ 

TABLE 1. Volume transport above 600 m (10 m^/sec) 
and distance AX west of the 1000 fm curve. 
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C. CURRENTS AT A DEPTH OF 200 METERS 

At a depth of 200 meters (Figures 17-22) current speed 
and direction were remarkably similar to those found at the 
surface on all six cruises with few notable exceptions. 
During the winter season (Cruise III) the offshore surface 
anticyclonic eddy reversed its direction to cyclonic at 
200 m. In addition, the eddy was found in the middle of a 
large band of poleward flow Instead of between separate 
fingers on northerly and southerly flow as was the case on 
the surface. The only other notable exception occurred in 
the early summer cruise (Cruise V) where the northernmost 
eddy reversed directions between the surface and 200 m. 

It should be noted that these apparent eddy reversals may 
be a distortion of the geostrophic flow due to the influence 
of internal waves . 

D. TEMPERATURE AND SALINITY AT THE SURFACE 

Figures 23, 2^1, and 25 show that tongue-like distribu- 
tions of temperature and salinity at the sea surface were 
encountered during all seasons of the year (Cruises III, V, 
and VI). During the winter season (Cruise III) a tongue of 
cold (less than 13°C) relatively fresh water was found 
stretching to the southeast between ^l0-60 nm offshore. A 
tongue of water with salinity less than 31 Voo was also 
found near the mouth, and to the north of, Monterey Bay. 

A comparison v;ith the surface flow (Fig. 13) shows that the 
tongue of cold relatively fresh water coincides with the 
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Fig. 24. Cruise V. Temperature (°C) and salinity (%o) at the sea surface. 
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Fig. 25. Cruise VI. Temperature ' (°C) and salinity (%©) at the sea surface. 



offshore band of current flowing to the southeast. The 
nearshore tongue of low salinity water coincides with the 
nearshore cyclonic eddy which perhaps indicates the eddy 
may not be closed as conjectured, but may be a finger of 
the current from the north. 

During the early summer cruise (Fig. 2^) a band of low 
surface temperatures (less than 11°) stretched seaward from 
the coast to a distance of more than 15 miles. A tongue of 
water with salinities higher than 33. 8% o was found in the 
same region. These low surface temperatures and high 
surface salinities are most likely due to upwelllng in the 
region. A comparison with the surface flow (Fig. 15) shows 
that the tongues of higher temperature coincide with the 
offshore bands of poleward flow. 

Conditions found during late summer (Cruise VI) were 
somex^hat similar to those found during early summer in that 
a tongue of relatively cold, high salinity water was encoun- 
tered near the coast off Monterey Bay. As will be shown in 
later sections, this was apparently due to intermittent 
•upwelllng and indicates that the upwelllng season had not 
completely ended. Hov/ever, it is apparent from the consid- 
erably higher offshore temperatures that a transition to the 
oceanic season was taking place. A comparison v;ith the 
surface currents (Figure l6) shows that the two eddies 
appear to be associated with areas v/here extreme bending of 
the isotherms occur, and v/ith areas of relatively low 
salinity . 
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E. TEMPERATURE AND SALINITY AT 200 METERS 

At a depth of 200 m (Figures 26, 2? and 28) the tongue- 
like projections of isotherms and isohalines, similar to 
those found at the surface, were noted during the early and 
late summer cruises. However, during the winter season 
(Cruise III) the isolines of temperature and salinity were 
more nearly parallel to the coast . The band of salinity 
minimum that bisected the region shows, upon comparison 
with the 200 m flow (Figure 19), that it conformed to the 
narrow finger of equatorward current. The offshore cyclonic 
eddy was nestled in the area of divergence of the 8.7®C 
isotherms, whereas the nearshore eddy coincided with a 
region of minimum temperature. It is believed that this 
small area, as well as the region of relatively cold, high 
salinity water near Pt. Sur were perhaps due to intermittent 
upwelling. This conclusion will be substantiated in later 
sections . 

During early summer (Pig. 27) the extreme bending of 
the 7.6®C and 7.8®C isotherms intersected with the tongues 
of highest salinity and coincided with relatively high 
speed poleward flov; nearest the coast (see Pig. 21). As 
will be shown in later sections, this was apparently a 
region of maximum intrusion of Equatorial Pacific v/ater. 

During late summer (Fig. 28) the lowest temperatures 
and salinities v^ere found farthest from shore. A comparison 
with the 200 m flov; (Fig. 22) shows that the offshore 
tongue-like intrusion of cold, relatively fresh v/ater 
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Fig. 26 . Cruise III. Temperature (°C) and salinity (%o) a-t 200 
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Fig. 27. Cruise V. Temperature (°C) and salinity i%o) at 200 
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Fig. 28. Cruise VI. Temperature, (°C) and salinity (%o) at 200 



coincides with the offshore southward current. As in the 
early summer months (Cruise V) the nearshore region off the 
northern edge of Monterey Bay was also characterized by 
relatively cold, high salinity water apparently due to 
upwelllng. Also of interest Is that In all three cruises 
the areas of upwelllng were the same areas where eddies 
were found. 

The agreement between the currents and patterns of 
temperature and salinity was remarkable In view of the 
probable non-stationarity and the time Interval between 
observations. The agreement may be interpreted as indicating 
that, even with non-stationary conditions, the changes were 
slow and modified only the details of the general picture, 
especially at locations away from the shelf. 

F. VERTICAL DISTRIBUTION OF TEMPERATURE AND SALINITY 

The variations in the vertical distribution of tempera- 
ture and salinity between the winter season (Cruise III), 
early summer (Cruise V) and late summer (Cruise VI) are 
depicted graphically in Figures 29-3^. A comparison of the 
temperature profiles (Figs. 29, 31 and 33) shows that the 
thermocline extended from the surface to nearly 100 m at 
the stations farthest from shore during late summer, and in 
contrast to the winter and early summer cruises, it inten- 
sified as it approached the coast, rising to a thin surface 
layer of only 25 m at the nearshore stations. During all 
seasons the surface temperature increased with distance 
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DISTANCE FROM SHORE (NM) 

Fig. 31. Cruise V (June-July 1971). Vertical cross section of temperature 
between the surface and 500 meters. The section runs WSW from Monterey Bay. 
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from shore, a phenomenon which Is typical of the region 
according to Reid (i 960 ). In addition, the seasonal vari- 
ation of the thermal structure confirms the annual cycle 
postulated for this region by Skogsberg (1936). 

Marked variability in salinity distributions between 
seasons can be seen in Figures 30, 32. and 3^. The most 
significant change between seasons can be described simply 
by looking at the 3^.0 %© isohallne . Note that in winter 
it remains constant at approximately 200 m, whereas in late 
summer it rises rapidly from a depth of 3 OO m at the farthest 
offshore station, to a depth of about 60 m approximately 
5 nm from shore. In early summer it, as well as the deeper 
isohallnes, rises significantly in an area about 25 nm 
offshore. As will be shown in a later section concerning 
distribution of water masses, this rise was apparently due 
to an intense intrusion of Equatorial Pacific water. 

The 26.4 isopycnal is- illustrated (dashed line) on 
Figures 30, 32 and 3^ to show the results of the above 
distributions. It is nearly , horizontal and constant at 
about 200 m in winter. However, in late summer it rises 
from a depth of more than 240 m offshore to less than 75 m 
at the nearshore station. Comparing Figures 25, 28 and 3^ 
confirms the previous speculation that late seasonal up- 
welling has occurred in the nearshore region during the late 
summer of 1970. It should be noted that the time frame 
during which the upwelling actually occurred is very diffi- 
cult to define. It can only be concluded that either 
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upwelllng was actively occurring during the time of the 
cruise, or the water was from a previous upwelllng event 
which had been horizontally advected Into the area. 

G. MIXING OF WATER MASSES BETWEEN 100 AND 1000 METERS 

As previously discussed in Section II, two water masses 
are typically found in the California region, a northern 
(Subarctic) water mass and a southern (Equatorial Pacific) 
water mass. These water masses are not characterized by a 
single point on a T-S diagram, but rather by T-S curves. 

In any consideration of the mixing of water masses charac- 
terized by T-S curves certain assumptions concerning the 
character of the mixing must be made (Sverdrup and Fleming, 
19^1). The assumptions made for this study were that the 
observed distribution of temperature and salinity in the 
area of transition is maintained by essentially horizontal 
advectlon and by lateral mixing of waters having the T-S 
relationships defined in Figure 3, and that mixing takes 
place primarily along surfaces of equal a^. On the basis 
of these assumptions a diagram similar to that used by 
Tlbby (19^1) was constructed (Figure 35). The area between 
the two T-S curves, defining the Subarctic and the Equatorial 
Pacific v;ater masses, has been divided into proportioned 
parts, and new T-S curves have been constructed which corre- 
spond to mixtures of Subarctic water containing 10 to 90 
per cent of Equatorial Pacific v;ater. By means of this 
diagram the percentage of Equatorial Pacific v;ater which 
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Pig. 35. Graph showing T-S curves defining Subarctic 
water and Equatorial Pacific water, and curves for 
various percentages of Equatorial Pacific water assuming 
mixing along surfaces of equal a^. 
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corresponds to any given temperature and salinity within 
the entire area examined can be charted and studied. It 
should be noted that the determination of percentage compo- 
sition by this means cannot be used for water above about 
100 m, since near the surface, non-conservative processes 
are prominent and mixing along surfaces will be masked 
(Tlbby, 19^1). Below a depth of about 1000 m the differences 
in T-S relationships of the two extreme water masses are so 
small that application of the method is uncertain. Figures 
36 , 37 and 38 , which are T-S diagrams typical of the winter 
season, early and late summer, respectively, point out these 
facts graphically. Note that on the diagrams symbols are 
placed at the 100, 200, 300, il00, 500, 6 OO, 8 OO and 1000 m 
levels. The two curves in Figure 38 illustrate the maximum 
and minimum T-S values found below 100 m during all three 
seasons. Only Station 10 on Cruise V (early summer) had 
values outside the range established by the two curves. 

Figures 39 » ^0 and hi are vertical sections showing the 
percentage of Pacific Equatorial water for the winter season, 
early and late summer, respectively. One feature of the 
distribution is common to all seasons: the percentage is 

higher tov^ards the bottom and towards the shore. This 
feature was found to generally hold true for all six cruises 
and agrees with the results obtained by Tibby (19^1). As 
can be seen in Figure h2, his results differ from those 
obtained in this study in that he found a smaller percentage 
of Equatorial Pacific v/ater in the area. Figure h2 shows 
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Fig. 36 . Cruise III (Feb. 1970). T-S curves from 
the surface to 1000 m for stations Q and R. 

Symbols indicate selected depths (see text). 
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Fig. 37. Cruise V (June-July 1971). T-S curves 
from the surface to 1000 m for Stations 1-8 and 
I-IO. Symbols indicate selected depths (see text). 
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Pig. 38 . Cruise VI (Aug. 1970). T-S curves from 
the surface to 1000 m for Stations H and 6. 
Symbols Indicate selected depths (see text). 
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Fig. 42. Section 3 off Monterey, California, showing 
percentage of Equatorial water (from Tibby, 19^1). 



76 



that he found a maximum of 50 per cent Equatorial Pacific 
water In the area, whereas greater than 60 per cent Equatorial 
Pacific water was found on all six cruises examined In this 
study. 

Figure 40 shows a significant southern water intrusion 
of considerable vertical extent near Station 10. When 
compared to Figures 21, 2? and 32, It confirms the earlier 
assertion that the distribution of properties In this area 
were due to the Intensity of the Equatorial Pacific Intrusion. 

Figures 43 and 44, in addition to Figure 4l, show the 
latitudinal variability of percentage composition of Equa- 
torial Pacific water found In late summer (Cruise VI) . 

Figure 43 shows the northernmost line while Figure 44 shows 
a transect about' 60 nm to the south. A comparison of these 
two figures with the flow at the surface and at 200 m (see 
Figs. l6 and 22) shows very good agreement. The offshore 
stations having the low percentage of Equatorial Pacific 
water are being fed by the California Current carrying 
Subarctic water. The nearshore stations lie in the path of 
the Davidson Current, thus giving rise to a greater percentage 
of Equatorial Pacific water. 
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VI. SUMMARY AND CONCLUSIONS 



The oceanic region off the Central California coast 
between Cape Mendocino and Pt. Conception has received 
relatively little attention In the literature since the 
early 19^0' s. A knowledge of the physical oceanography 
of this region Is of vital Importance to commercial enter- 
prises In the area and to new projects such as the construc- 
tion of a nuclear power plant near Davenport, California. 

This study has attempted to describe the seasonal vari- 
ations In currents, temperature and salinity, and water 
masses In a portion of this unexplored region. The data 
discussed were obtained on six cruises covering a span of 
two years from 1970 to 1972. 

A. CURRENTS 

Computations based on density distribution show that 
the surface currents are subject to marked seasonal varia- 
bility. The area studied Is generally characterized by 
small scale eddies, especially during the upwelllng season, 
and by the Interlacing of separate bands of poleward and 
equatorward flov;. Assuming conservation of vortlclty Is 
valid In this region. It can be concluded that the bottom 
topography Influences the direction of flow Inside the 1000 
fathom curve . 

At a depth of 200 m the currents are very similar to 
those found at the surface. 



80 







4 







► 



a 







B. TEMPERATURE AND SALINITY 

Skogsberg's ( 1936 ) annual cycle, composed of three 
distinct hydrographic seasons, was generally adhered to 
during the two year span that this study covered. The 
results confirm the existence of a winter and upwelllng 
season, Although the existence of a distinct oceanic 
season could not be confirmed due to a paucity of data, 
the results of the late summer cruise (Cruise VI, August 
1970) clearly show a transition towards the conditions 
characteristic of an oceanic season. 

C. MIXING OF WATER MASSES 

The area studied lies virtually In the middle of the 
transition zone between Subarctic and Equatorial Pacific 
waters. As a consequence, the mixing of these two extreme 
water masses Is nearly complete. Thus only subtle changes 
from season to season are noted, with the highest percentage 
of Equatorial Pacific water being found. In general, during 
winter and early summer. The results of this study confirm 
Tlbby's (19^1) results In that the percentage composition 
of southern water Increases v;lth depth and towards shore. 
Greater than 60 per cent Equatorial Pacific viater Is found 
throughout the year below 8OO m. 
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APPENDIX A 



"HYDRO" 

A Computer Program for Processing 
Hydrographic Data 



A . INTRODUCTION 

This discussion is presented in order to familiarize 
the reader with a computer program that will Interpolate 
salinity and temperature data for "standard depths" and 
calculate density (sigma-t), specific volume, specific 
volume anomaly, sound velocity, dynamic height, dynamic 
height anomaly, and geostrophlc currents and transports. 

The theoretical development of the appropriate equations 
may be found in Sverdrup, Johnson, and Fleming, The Oceans , 
or in Neumann and Pierson, Principles of Physical Ocean- 
ography . The application of these equations is embodied 
in the text of the program. 

B. PREPARATION OF THE DATA DECK 

The data deck consists of three main parts: 

(a) number of geostrophlc current calculations to be 
performed 

(b) designation of pairs of stations between which 
geostrophlc currents and transports are to be 
computed 

(c) station data. 

These parts are discussed individually in more detail 
assuming the user has a fundamental knowledge of computer 
programming in the Fortran language . 
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(a) Number of Geostrophlc Calculations (termed NGC) 
refers to the number of pairs of stations between which 
user wishes to have geostrophlc calculations made. This 
two digit Integer Is punched In columns 3 and 4 on the 
first data card as Indicated below. If no geostrophlc 
calculations are desired, the first data card should have 
zeros (or blanks) In columns 3 and 4; and the user should 
skip to part (c), preparation of station data. 



Column 


Format 


Item 


1-2 


2X 


blank 


3-4 


12 


NGC (maximum value Is 48) 



(b) If NGC was not set equal to zero, the next card(s) 
must specify NGC pairs of stations and the accompanying 
reference levels for the geostrophlc currents. 

At this point the user will find It helpful to plot his 
hydrographic stations on a chart of some sort. This chart 
should be sufficient to Illustrate spatial relationships 
between stations . The user should also determine from the 
chart the sequence In which the stations will appear In 
the data deck. 

It Is helpful to think of the hydrographic data as 
being stored In a two dimensional array. The station ap- 
pearing first In the data deck Is stored In column 1, the 
second station In column 2, and so on. Within each column 
the temperature, salinity, and oxygen data are stored 
according to depth. The specification of a pair of stations. 
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station A and station B, is made using these column index 
numbers. These are the sequence numbers determined from 
the station plot. 

The signs of the calculated relative and absolute veloc- 
ities should be interpreted according to the following con- 
vention. If station A is chosen to be to the right of 
station B, and the dynamic height at station A is greater 
than the dynamic height at station B, then the relative 
velocity is positive directed away from the observer. That 
is, the quantity (v^ - V2) is greater than zero, where v^ 
refers to a surface current and V2 refers to a deeper current 
as shown in the figure. 



« B sea surface A . 




If > Dg, 

then “ ^2) > 0 . 



(Note: If it can be assumed that the current is dominated 

by a north-south component, then station A should be chosen 
so that it is east of station B. If it can be assumed that 
the current is dominated by an east-west component, then 
station A should be chosen so that it is south of station 
B. This technique v/ill assure that the calculated currents 
possess signs according to the generally accepted convention: 
V positive northward and u positive eastward.) 



The reference level must be one of the standard depths 
defined in the program: 0, 10, 20, 30, 50, 75 > 100, 125, 

150, 200, 250, 300, il00, 500, 600, 700, 800, 1000, 1200, 
1^00, l600, 2000, 2500, or 3000 meters. (These depths can 
be redefined if a particular situation requires it, simply 
substitute the appropriate card.) If the user specifies a 
reference level deeper than the data for the two stations, 
the program will automatically use the deepest level common 
to both stations. 

The format for the second part of the data deck is as 
follows : 



Column Format Item 



1 


IX 


blank 


2-3 


12 


station A 


H-5 


12 


station B 


6-10 


P5.0 


reference level 


11 


IX 


blank 


12-13 


12 


station A 


14-15 


12 


station B 


16-20 


P5.0 


reference level 



. . . and so on until NGC pairs of stations have been 
specified. Recall that data cards are read from column 1 
through 80. 

(c) The data for each station is then presented as 
follows : 

i. A heading card that indicates the number of 

measurements (NOV) of temperature, salinity, and 
oxygen taken at the station, the station name 
(for the user’s convenience only), latitude. 
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longitude, and the date/time group. It would be 
helpful to suffix the station name with the 
sequence Index number If It Is known at the time 
the data deck Is prepared. 

11. NOV data cards that list temperature, salinity, 
and oxygen values as a function of Increasing 
observed depth. Depth values here are specified 
as positive quantities. Oxygen values are not 
used In any program calculations and may be left 
blank on the data cards. 



Heading card: 






Column 


Format 


Item 


1-H 




NOV (maximum value Is 50) 


5-16 


3Ai| 


station name 


17-19 


F3.0 


latitude (degrees N) 


20-2^ 


F5.1 


latitude (minutes to nearest 






tenth) 


25-28 


Fil .0 


longitude (degrees) 


29-33 


F5.1 


longitude (minutes to nearest 






tenth) 


3iJ-iJ5 


3Ail 


date/time group 



Station data card: 



Column Format Item 



1-8 


F8.1 


observed depth 


9 


Al 


(see below) 


10 


IX 


blank 


11-19 


F9.2 


temperature (°C) 


20 


Al 


(see below) 


21-30 


F10.3 


salinity (°/oo) 


31 


Al 


(see below) 


32-39 


F8.2 


oxygen (ml/1) 


^0 


Al 


(see below) 




8x 


blank 


^49-80 


44A8 


comments that user wishes to 



Include; these comments will 
appear alongside station data 
In the printed output 
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The alpha space immediately behind the values of depth, 
temperature, salinity and oxygen content may be used to 
indicate the supposed validity of that value . For instance , 
if the value is questionable, a question mark (?) may be 
used; if the value is adjusted relative to the other data, 
an asterisk (*) may be used; and if one wishes to indicate 
no value, the letter X may be used. This last notation is 
particularly useful in the oxygen column where some values 
might be missing and the normal "0.0" output would be 
misleading. 

(d) After all station data have been assembled, the 
final card in the data deck must be a blank. 

C. CALLING SEQUENCE FOR HYDRO 

Program HYDRO is currently stored as an object deck 
on a disk at the NPS computer facility. This program is 
in machine language in order to minimize compile time. 

The source language version of the program is presented 
at the end of this instruction. 

If an extenuating case exists, the source deck for 
HYDRO may be obtained from the Oceanography Department. 
However, for most cases it is advantageous for the user to 
access the library version as follows: 
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# 




//Job card 

//J0BLIB DD DSNAME=F2853.0CEAN0,DISP=SHR, 

// V0LUME=SER=CELOO2,UNIT=2321 

//G0 EXEC PGM=HYDR0,REGI0N=1OOK 

//PT 06 P 001 DD SYS0UT=A,SPACE=(CYL,(1,1)), 
// DCB=(RECPM=PBA,LRECL=133,BLKSIZE=3325) 

//PT05F001 DD * 

(data deck as described above) 

/* (standard orange card) 
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REAL *8 ITL ( 12 ) , INFO( 50,4) 

DIMENSION 10(50) , IT (50), IS (50), 10(50) 

DIMENSION 0(50) ,T(50) ,S(50) ,02 (50) ,SNDV(50) 

DIMENSION S0(24) ,ST(24) ,SS( 24) ,SGT(24) , SV( 24) , SVA( 24) 

DI MENS ION SGP( 24 ) , DH( 24) , b0H(24) ,00( 24) , SL EV( 48) , BSVA( 24) 

DIMENSION NPA(43) ,MPG(43) ,N0(43) ,NSTA(43,3 ) , ALT(40 ) ,AL:1(43) 
DIMENSION ALN (48 ) , ANN (43) , I DATE ( 48, 3) ,DHT( 48 ,24) ,A0H(24) 

8 FORMAT (1H1,« STATION, ' ,3A4, • LATITUDE ',I2,F5.1,'N LONGITUDE 

113, F5.1,'W DATE- • ,3A4/5X,' T0,',3A4,' LATITUDE *,I2,F5.1, 

,2'N LONGITUDE • , 1.3, F5e 1 , ' W DATE ',3A4/./) 

9 FORMAT (8(1X,2I2,F5.0) ) 

10 FORMAT ( IHl, ' STATION *,3A4,' LATITUDE = ',I2,F5.1,'N LONGITUDE 

114, F5.1,‘W DATE ' , 3A4// ) 

11 FORMAT (10X,'=- INDICATES ADJUSTED VALUE') 

12 FORMAT (lOX, '? INDICATES QUESTIONABLE VALUE') 

99 FORMAT (10X,'X INDICATES NO VALUE') 

13 FORMAT ( I4,5A4,F3.0,F5.1 ,F4.0,F5.1 ,3A4) 

15 FORMAT (F8.1,A1, 1X,F9.2,A1,F10.3,A1,F9.3,A1,7X,4A8) 

16 FORMAT (lOX, 'DEPTH TEMPERATURE SALINITY SIGMA-T OXYGEN'//) 

17 FORMAT ( 20X, ' OBSERVED VALUES'//) 

18 FORMAT (//20X, ' INTERPOLATED VALUES'//) 

19 FORMAT (lOX, 'DEPTH TEMPERATURE SALINITY SIGMA-T SNDVEL SPE 

1C VOL SPEC V ANOM MEAN SVA DELTA 0 DYNAMIC HEIGHT'//) 

20 FORMAT (1 OX, F5.0, F10.2, F12. 3 , F9 . 3 , F 12 . 2, F 10. 4, F 12. 6 , 25X, F 1 0. 5/8 0X , 
12F12.6) 

21 FORMAT (10X,F6.1,A1,F8.2,A1,F11.3,A1,F8.3,F9.3,A1,4X,4A3/) 

-DATA SD/0. ,10. , 20. ,30. ,50.., 7 5..., 100.., 125. , 150. ,200.. , 250. , 300. ,4 00. , 
1500. ,6 00. ,700. ,800. ,1000. ,1200. , 1400 ., 1600 ., 2000 . , 2500. ,3 000./ 

READ THE NUMBER OF 6E0STR0PHIC CURRENTS C TRANSPORTS 
TO BE CALCULATED 

READ (5,13) NGC 
IF(NGC.EQ.O) GO TO 410 

READ (5,9) (NPA( I ) ,NPB( I) ,SLEV( I) , I=1,NGC) 

410 NPA(NGC+1)=0 
DO 41 L=l,48 

READ HEADING CARO, CHECK FOR END OF DATA, THEN 
READ NOV DATA CARDS. 

READ (5,13) NOV, (NSTA( L,K) , K=1 , 3 ) , ALT( L ) , A L M ( L ) , ALN ( L ) , ANM ( L ) , 
1(I0ATE(L,K),K=1,3) 

IF (NCV) 32,32,24 
24 DO 25 1=1, NOV 

READ (5,15) 0( I ) , I0( I ) ,T(I ) , IT( I ),S( I ) ,IS( I ) ,02(1 ) , 10(1 ) , 

1(1NF0( I ,JJ , J = 1 ,4) , 
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SGTSVA IS SUBROUTINE TO COMPUTE SIGMA-T, SPECIFIC VOLUME 
AND SPECIFIC VOLUME ANOMALY. 

25 CALL SGTSVA ( T ( I ) , S ( I ) , D ( I ) , SGP ( I ) , S VNO, S VNO ) 

LGTP IS SUBROUTINE TO COMPUTE INTERPOLATED VALUES 

CALL LGTP(NOV,D,T,SD, ST,NA) 

CALL LGTP(NOV,D,S,SD,SS,NB) 

..,NO(L)=NA 
DO 27 1=1,NA 

SNDVEL IS SUBROUTINE TO COMPUTE SOUND VELOCITY 

CALL SNDVEL (ST ( I ) ,SS ( I ) ,S D( I ) , SNDV( I ) ) 

CALL SGTSVA ( S T( I ) , SS ( I ) ,S D ( I ) , SGT ( I ) , S V ( I ) , S V A ( I ) ) 

27 CONTINUE 
NLT=ALT(L) 

NLN=ALN( L) 

WRITE (6,10) (NSTA(L,K),K=1,3),NLT,ALM(L) ,NLN,ANM(L) , 
1(IDATE(L,K),K=1,3) 

WRITE( 6, 17) 

WRITE (6,16) 

DO 29 I =1 ,NCV 

29 WRITE (6,21) D ( I ) , I D ( I ) , T( I ) , I T ( I ) , S ( I ) , I S ( I ) , SGP ( I ) , 02 ( I ) , 1 0( I) , 
1(INFC(I,J),J=1,4) 

WRITE (6,11) 

WRITE( 6, 12) 

WRITE (6, 99) 

WRITE (6,18) • 

WRITE! 6, 19) 

NA=NA-1 

DH(L)=0. 

DO 30 I=-'1,NA 

BSVA( I ) = (SVA( I H-SVA( I + l) )*.5 
DD( I ) =BSVA(I )=^ (SD( I + l )-SD( I ) ) 

30 DH( I+l )=DH( I )+D0( I ) 

DO 31 1=1 ,NA 
DHT(L, I ) =DH( I ) 

31 WRITE (6,20) SD ( I ) , ST ( I ) , SS ( I ) , SGT ( I ) , SNDV ( I ) , S V ( I ) , SVA ( I ) , DH( I ) , 
1BSVA( I ) ,DD( I ) 

I=NA+1 

DHT(L, I ) = DH( I ) 

WRITE (6,20) SD( I ) ,ST( I ) ,SS ( I ) , SGT ( I ) , SNDV( I ) , SV( I ) , SVA( I ) ,DH( I ) 
41 CONTINUE 

IF(NGC.EQ.O) GO TO 33 

32 DO 42 L=1 ,48 

IF (NPA(L) .EQ. 0) GO TO 33 
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I 






BASE=SLEV(L) 

N1=NPA(L) 

N2=i\'PB(L) 

NU1=N0( M) 

NU2=N0(N2) 

DO 43 l=l,i\Ul 

43 AOH( I ) =DriT(Nl , I ) 

DO 44 1=1, KU2 

44 BDH{ I ) = CriT (N2, I )• 

NLT=ALT(.\1) 

NLN=ALN(N1) 

MLT=ALT(.N2) 

KLN=AL,\(.N2) 

WRITE (6,8) (NSTA{N1,K) ,K=1,3) ,NLT,ALfUNl) ,NLN,ANM(N1) , 

1 (IDATE{.vi ,K) ,K=1,3 ) , (NSTA(iM2,K ) ,K=1, 3) , MLT , ALM ( N2 ) , flLN , 
2ANM(iM2) ,(IDATE(N2,KJ ,K=1,3) 

ALAT=ALT (M1)+ALM(M1 )/60. 

ALON=ALi\(Nl l + ANrUi-a )/60. 

BLAT=ALT( N2 ) +A LM ( N2 ) /60. 

B L 0 N = A L N ( !nI 2 ) + A N M { N 2 ) / 6 0 . 

DSTSTA IS SUBROUTINE TO COMPUTE DISTANCE BETWEEN STATIONS 

CALL DSTSTA { AL AT , ALON , BLAT , BLON, X2, DI ST ) 

GEOCUR IS SUBROUTINE TO COMPUTE GEOSTROPHIC CURRENTS 
AND -TRANSPORTS 

42 CALL GEOCUR { NUl , ADH, NU2 , BDH , SD, BASE , X2 , NNN , D I ST) 

33 STOP 
END 
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SUBROUT IME LGTP(N,D, V, SD,CV,NN) 

01 MENS ION D(N) ,V (N ), CV(24) ,SD( 24) 

111 DO 183 J=l,24 

112 DO 186 1=1, N 

IF(SD( J)-D(N) ) 113, 115, 190 

115 CV(J)=V(N) 

GO TO 191 

113 IF(SO( J)-0(I))114,114,116 

114 CV( J)=V( I ) 

•GO TO 138 

116 IF(SD( J)-D(H-l) }120,118,186 
118 CV( J)=V( I+l) 

GO TO 133 

120 IF(I-l) 132,132,126 

126 XA=(SD(J)-D( I) )^=(SD(J)-D(I + 1))«V(I-1)/ 
1( (D( I-13-D( I ) )=MD( I-l )-D( I + l) ) ) 

XB=( SD(J)-D( I-l) )*(SD( J)-D( I+l) )^V(I)/ 
1( ( 0( I )-D( I-l) )*(D( I )-D( I+l) ) ) 
XC=(SD(J)-D( I-l ) )=i^(SD( J)-D( I ) )*V (I + l)/ 
l((D(I+l)-0( I-l) )*(D(1+1)-D( 1) ) ) 
ANSU=XA+XB+XC 

132 IF ( (I+2)-N) 133,133, 134 

133 YA=( SD(J)-D( I+l) )*(SD(J)-0(I+2) )-V(I)/ 
1( (D( I )-D( I+l) ) *( D( I )-0( 1+2) ) ) 

YB=( SD( J)-D( I ) )=?(SD( J)-D( 1+2) )*V(I+1)/ 
1( (D( I+l)-D( I ) )*(D( I+l)-D( 1+2) ) ) 

• YC=(SD(J)-0( I ) )*(SD(J.)-D( I + l) )-V(I+2)/ 
1( (D( l+2)-D( I ) )-(D( I+2)-D( 1+1) ) ) 
ANSD=YA+YB+YC 

134 ZA=(SD(J)-D( I+l ) )*V(I )/(D( I )-D( I+l) ) 
ZB=( S0( J)-D( I ) )=^'V( I+l )/(D( I+l ) -D(I ) ) 
ANSL=ZA+ZB 

IF(I-1 )136,136, 138 
136 CV( J)= (ANSD+AMSD/2. 

DLL=( ANSD+ANSL+ANSL)/3. 

GO TO 188 

138 IF( ( I+2)-N) 140, 140,142 
140 CV( J )= (ANSU+ANSD+ANSD/3. 

UD=( ANSU+ANSD)/2. 

GO TO 188 

142 CV ( J ) = ( ANSU + ANSD/2. 

ULL=( ANSU+AMSL+ANSL)/3 . 

GO TO 188 
186 CONTINUE 
188 CONTINUE 

190 J=J-1 

191 NN=J 
RETURN 
END 
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SUBROUTINE SNDVEL (T,StD,SVU 
IF(D.LE.100.) GO TO 10 
IF(D.LE.200. ) GO TO 20 
IF(D.LE.400. ) GO TO 30 
IF(D.LE.700.) GO TO 40 
IF(D.LE.1500. )G0 TO 50 
IF(D.L£.2000. )G0 TO 60 
P=D=:'-{ .10318-M0-2000.)*2.6E-7) 

GO TO 5 

- 1.0 P=D*( . 10245+D’!=3.0E-7) 

GO TO 5 

20 P=D-(.10248+(0-100.)*7.0E-7) 

GO TO 5 

30 P=D*( . 10255+{D'-200.)*6.0E-7) 

GO TO 5 

40 P=D*{ . i0267+(D-400. )=;=4.33E-7) 

GO TO 5 

50 P=D*(. 10280+(D-700. j~3.0E-7) 

GO TO 5 

60 P=D*{ .10304+(D-1500. )*2.7E-7) 

5 CONTINUE 

VT=T=:={4. 5721-T*( . 044532-T* ( 2. 6045E-4^•T* ( 7. 9851 E-6 ) ) ) ) 
VP=P*( .160272 + P*(1 .0268E-5+P*(3.5216E-9-P=i={3.3603E-12) ) ) ) 
VS=( S-35. ) *{ 1.39799+( S-3 5 . ) *1 . 69202 E-3 ) 

A1=T^M-. 011244+3. 1530E-8^-p + T^={7.7711E-7+l. 5790E-9*P) ) 

A2 = P-(7.7016E-5-1.2943E-7='=P) 

A3=--P*T-M-1.S607E-4+T^‘= ( 7. 4312E-6+P" {4.5233E-8 ) ) ) 

A4=P*P*(T*{-2.5294E-7+T-1. 8563E-9)+P"T^M-1.9646E-i0) J 

SVL=1449. 14+VT + VP+VS+(S-35 . ) =? { A1+A2 ) +A3 +A4 

RETURN 

END 
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SUBROUTINE SGTSVA ( T , S , 0, SGT, SV, SVA ) 

ST=-( ( (T-3.98)^=’^2) /503.57)*( (T+283 . ) / ( T+67 .26 ) ) 
CL=(S-.030)/l. 305 

S0=-.069+1.4708*CL-.00157«CL**2+3,98E-5*CL**3 

ALTERNATE METHOD OF COMPUTING SIGMA-ZERO 
S0= -0.093 + 0. 8149*S-.000482>«^S'+=!=2 + 6. 8E-6*S 

AT=T*{ 4.7367-. 0981 3 5*T+. CO 1 0843*T**2 ) * 1 . E-3 
•BT=T*( 13 .030-. 5i64=!=T+..0.1.66 7=i--T-’’-2 )*l.E-6 
SGT=ST+(Sa+. 1324)-(l.-AT+oT*(S0-.1324) ) 
AFST=l./( l. + SGT*l.E-3) 

A=D=:‘AFST*l.E-9 
B=4886. /(l.+l. S3E-5*D) 

C=227. + 2S.33-T-.551=!^T**2 + . 004« 

E=D«l.E-4 
G=(S0-28. ) /lO. 

H=147. 3-2.72*T+.04*T«*2 
U=10 5 . 5 +9 . 5 -T- . 1 5 8 *T *=^2 
V=l. 5*0— ^2 — T— 1 .E — 8 
W=32.4-.87^T+.02-T-*2 
X=4.5-.1*T 
Y=l. 8- . 06=i'T 

SV=AFST-A=!^(B-C + E*U-V-G=?=( H-E*W) +G*=<'2=<'{ X-E=!=Y) ) 
AZ=. 972843 
YA=-227.+. 01055*0 
YB=. 0126-'-( 147. 3-^.00324*0) 
AP=AZ-D*AZ*(B+YA-YB)*l.E-9 
SVA=SV-AP 
RETURN 
END 
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SUBROUTINE D5T ST A( S AT I , ONG I , S AT 1 1, ONGI I , X2 , 0 1ST ) 

IMPLICIT REAL^4 (K) 

‘ REAL-8 A,E 

DATA A/111132. 09/ , B/566 . 05 / , C/ 1 . 20/ , 0/ .002/ 

DATA E/ 111415. 13/, F/94. 55/,G/.012/ 

10 FORMAT (lOX.'HEAN LATITUDE = ' , F6. 2/ 15X, 'D I STANCE = SF6.2, 
1* KILOMETERS’/) 

C0N=2«3. 1416/360 
AATI=SAT1*C0N 
. ..AATI I = SATI I*CCM 

$MERI = A-S-COS( 2-AATI ) +C=>COS( 4*AATI )-0*C0S( 6=«^AATI ) 
PARI=E*CGS( AAT I ) -F==^COS ( 3-A AT I ) +G^^C0S ( 5*AAT I ) 

$KERI I=A-B-CGS( 2^ AATI I ) +C-COS ( 4- AAT 1 1 ) -0=:=COS (6#AAT 1 1 ) 

PARI I = E«CCS( AATI I ) -F-COS ( 3*AAT 1 1 ) + G*COS( 5*AATI 1) 

ALLAT=($M£RI + 1..'1ERI I)/2 

ALLON=( PARI+PARI I ) /2 

DLAT=SATI-SAT I I 

DLON=ONGI^CNGI I 

KLAT=DLAT^ALLAT/1000 

KLONG=DLON'':=ALLON/1000 

KDIX = SQRT(KLAT^*2 + KLCNG=S'=i'2 ) 

DIST=KOIX 
W2=1.453E-4 
PSI=(SATI+SATI I )*0.5 
PS (2.-3. 14159/360. ) *PSI 
. SPSI = S If'KPSJ ) 

.IF-( SP-SI. LT.0.1 ) SPSI = 0.1 
X2=l./(W2-SPSr*KDIX) 

WRITE(6,10) PSI,KDIX 

RETURN 

END 
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SUBROUTINE GEOCUR ( NA , ADH , N5 , BDH , SD , BAS E, X2 ,NNN, DIST ) 

DIMENS lOfJ AOH(NA) , 6DH( N3) , SD( 24) ,RVEL ( 24) , VEL( 24) ,AMB(24) , AVT(24 ) 

10 FORMAT U3X,' DEPTH OYN HT DYN HT DIFF HT REL VEL AB S 

. IVEL ABS V0L'/14X,' M. STA A STA B A-B CM/SEC 

2 CM/SEC TRANSPORT -•/) 

11 FORMAT (13X,F5.0,2X,3( F9.5, IX) ,2 ( F8 .2 , 2X ) / 72X , F9 . 5 ) 

12 FORMAT ( • LEVEL OF NO MOTION MUST BE EQUAL TO A STANDARD DEPT 

IH * ) 

14 FORMATC • ,10X TOTAL VOLUME TRANSPORT IS COMPUTED BY SUMMING INCR 
lEMENTAL TRANSPORTS ABOVE LEVEL OF NO HOT I ON : • //5X , • TOTAL TRANSPORT 
2 PERPENDICULAR TO THE PLANE OF THE STATIONS IS ',F7.5,' SVERORUPS 
3RELATIVE TO ',F5.0,» METERS') 

15 FORMAT!//' * VALUES IN THIS COLUMN REPRESENT TRANSPORTS IN LAYER 
1 INCREMENTS'// ) 

IF(NA.LE.NB) GO TO 51 

N=NB 

GO TO 52 

51 N=NA 

52 DO 53 1=1, N 

AMB ( I ) =ADH( I )-EOH( I ) 

53 RVEL( I )=-AHB( I )*X2 
DO 54 1=1,24 

IF (BASE. EQ.SD! I ) )G0 TO 55 

54 CONTINUE 
WRITE (6,12) 

GO TO 70 

55 NM=I 

IF(NM.GT.N) NM^N 
BASE=SD(NM) 

DO 56 1=1, N 

56 VEL( I )=RVEL( I )-RVEL(NH) 

DO 553 I=2,N 

J=I-1 

AVEL=(VEL ( I)+VEL(J ) )="0 .005 
553 AVT( J)=AVEL-OIST=?(SD(I )-SD( J) )'-L.0E-03 
NM=NM-1 
VT=0. 

DO 57 I=1,NM 

57 VT=VT+AVT( I ) 

IF STATION A IS TO THE RIGHT OF STATION 3 AS AN OBSERVER LOOKS AT 
THE PLANE OF THE STATIONS, A POSITIVE CURRENT FLOWS AWAY FROM THE 
OBSERVER. 

58 WRITE (6,10) 

N=N7l 

DO 60 1=1, N 

60 WRITE (6,11) SD( I ) , ADH( I ) , BDH( I ) ,AMB( I ) ,RVEL( I ) , VEL( I ) ,AVT(I ) 
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N=N+1 
I = N 

WRITE (6,11) S0( I) ,AOH( I ) , BDH( I ) ,AMB( I )t RVEL( I ),VEL ( I ) 
NNN=N 

WRITE(6,15) 

WR1TE(6,14J VT,8ASE 
70 RETURN 
END 
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